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Abstract. First, we show that the gamma-ray bursts with dark optical afterglows (DOAs) cannot be explained 
by a failure to image deeply enough quickly enough, and argue that circumburst extinction is the most likely 
solution. If so, many DOAs will be "revived" with rapid follow up and NIR searches in the HETE-2 and Swift 
eras. Next, we consider the effects of dust sublimation and fragmentation, and show that DOAs occur in clouds 

of size R £ lOL^ 2 pc and mass M J> 3 x 10 5 L49 M©, where L is the luminosity of the optical flash. Stability 
considerations show that such clouds cannot be diffuse, but must be molecular. Consequently, we compute the 
expected column density distribution of bursts that occur in Galactic-like molecular clouds, and show that the 
column density measurements from X-ray spectra of afterglows, DOAs and otherwise, satisfy this expectation in 
the source frame. 



EVIDENCE FOR CIRCUMBURST 
EXTINCTION OF BURSTS WITH DARK 
OPTICAL AFTERGLOWS 

Optical afterglows have been detected for about 1/3 
of the rapidly-, well-localized gamma-ray bursts (e.g., 
[1,2]). This data-rich subsample of the rapidly-, well- 
localized bursts has naturally been the focus of the vast 
majority of the field's attention and resources over the 
past five years, but in the end, it is a biased sample. The 
nature of the dark optical afterglows (DOAs) has only re- 
cently become a subject of greater interest, and as might 
be expected given that these are by definition data-poor 
events, contradictory initial findings: [1] find that ;> 3/4 
of the limiting magnitudes for the DOAs are brighter than 
the detected, but faint, optical afterglow of GRB 000630, 
and conclude that the DOAs are consistent with a failure 
to image deeply enough quickly enough. However, [2], 
using a smaller sample of DOAs, find that the distribution 
of limiting magnitudes for the DOAs, even if treated as 
detections, is significantly fainter than the distribution of 
magnitudes of the detected optical afterglows, and con- 
clude that the DOAs cannot be explained by a failure to 
image deeply enough quickly enough. 

In this section, we first resolve this issue of whether 
the DOAs can be explained by a failure to image deeply 



enough quickly enough by applying Bayesian inference, 
a statistical formalism in which limits can be treated 
as limits, instead of detections (e.g., [3,4]). Next, we 
argue that extinction by circumburst 1 dust is the most 
likely explanation for the vast majority of the DOAs. 
Finally, we consider the effects of dust sublimation and 
fragmentation [7,8,9], and place constraints on the sizes 
and masses of the circumburst clouds of DOAs. 

Evidence that DOAs are Fainter than 
Detected Optical Afterglows 

The samples of [1] and [2] differ in the following 
ways: The sample of [1] is nearly complete - it con- 
tains limiting magnitudes for 95% of the DOAs preced- 
ing GRB 000630 - whereas the sample of [2] includes 
limiting magnitudes for BeppoS AX bursts only. Also, the 
sample of [2] relies less on GCN Reports, and more on 
published results. As BeppoSAX bursts appear to have 
been deeply imaged more often than bursts detected by 
other satellites, the different findings of these two papers 



By circumburst, we mean within the circumburst cloud, which is 
probably parsecs to tens of parsecs across ([5]; see also [6,7]). 
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FIGURE 1. R-band magnitude distribution of the detected 
optical afterglows (hashed histogram), and limiting magnitude 
distribution of the DOAs (unhashed histogram, added to the 
hashed histogram), scaled to 18 hours after the burst. The 
dotted curves are the best-fit models (see text) to these data. 
Clearly, most bursts have afterglows fainter than R ~ 24 mag 
18 hours after the burst, independent of the assumed shape of 
the brightness distribution. From [10]. 



tion of all of the bursts. The likelihood function is then 
given by 
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where N is the number of bursts in our combined sample, 
and 
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(e.g., [3,4]). We now consider a wide variety of two- 
parameter models for n(R): a Gaussian, a boxcar, an 
increasing power law, and a decreasing power law. We 
have fitted these models to the data in Figure 1, we 
plot the best-fit models in Figure 1 (dotted curves), and 
we list the best-fit parameter values and 68% credible 
intervals in Table 1 of [10]. Also in Table 1 of [10], 
we list the fractions /all of all bursts, and /doa of 
DOAs, that have afterglows fainter than R = 24 mag 18 
hours after the burst, the magnitude scaled to this time 
of the faintest detected optical afterglow in our sample. 
The results are fairly independent of the assumed shape 
of the brightness distribution: On average, we find that 
w 57±\\% of all bursts, or w 82+^% f DOAs, have 
afterglows that are fainter than R = 24 mag 1 8 hours after 
the burst. Consequently, the DOAs cannot be explained 
by a failure to image deeply enough quickly enough: 
As a whole, they are fainter than the detected optical 
afterglows. 



can be attributed to sample differences, and the practice 
of comparing limits to detections, which can turn such 
sample differences into sample biases, as we demonstrate 
below. 

But first, we combine the samples of [1] and [2], 
keeping only the bursts through GRB 000630 to maintain 
the completeness of the [1] sample. As is done in these 
papers, we scale the data to a common time: We use a 
temporal index of -1.4, the median temporal index in 
Table 1 of [2], and we scale the data to 18 hours after 
the burst, the median observation time in our combined 
sample. In the event of multiple limiting magnitudes for 
a single burst, we adopt the most constraining. We plot 
the combined sample in Figure 1 : The hashed histogram 
is a binning of the R-band magnitude distribution of 
the detected optical afterglows. The unhashed histogram, 
which is added to the hashed histogram, is a binning of 
the limiting R-band magnitude distribution of the DOAs. 

Although these distributions appear to be similar, since 
one is a distribution of limits, they can be consistent only 
if the vast majority of the DOAs could have been detected 
had they only been imaged a little more deeply. Being 
this unlucky, consistently, is of course very improbable. 
We quantify this improbability with Bayesian inference: 
Let n(R) be the normalized R-band magnitude distribu- 



Evidence for Circumburst Extinction and 
Against Alternative Explanations 

[11] find that the far-infrared luminosity-to-mass ratio 
for isolated and weakly interacting spiral galaxies is con- 
sistent with the average for Galactic molecular clouds, 
which suggests that most, if not all, of the star forma- 
tion in such galaxies occurs in molecular clouds. Since 
bursts likely occur where massive stars form (see, e.g., 
[12] for a review of some of the evidence), and since the 
optical depths through such clouds are typically many 
magnitudes (e.g., [5]), circumburst extinction is a likely 
explanation for the vast majority of the DOAs. 

[10] put circumburst extinction to a simple test: Con- 
sidering only (1) strongly collimated bursts that burn 
completely through the optical depth of their circumburst 
clouds (e.g., [7,8,9]), and (2) weakly collimated bursts 
that do not, [10] find that circumburst extinction can ex- 
plain snapshot optical vs. X-ray and optical vs. radio 
brightness distributions of all of the afterglows, DOAs 
and dark radio afterglows included. 

We now address a number of alternative explanations: 

Galactic extinction: [2] show that the distribution of 
Galactic column densities for the DOAs is consistent 
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FIGURE 2. Source-frame column density probability distri- 
butions Pi{N H ) for the four DOAs in Table 1 of [5]. The dotted 
curves mark bursts of unknown redshift, for which we adopt 
z = 1 (top panel), 2 (middle panel), and 3 (bottom panel). From 
[5]. 



with the distribution of Galactic column densities for the 
bursts with detected optical afterglows, ruling this out as 
an explanation for the vast majority of the DOAs. 

Host galaxy extinction unrelated to the circum- 
burst medium: If burst host galaxies are like our galaxy, 
this cannot explain the vast majority of the DOAs for the 
very same reason: Most lines of sight through our galaxy 
do not result in a great deal of extinction. 

However, the burst host galaxies might not be like our 
galaxy: [13] argue that the vast majority of the DOAs 
might be occurring in the nuclear regions of ultralumi- 
nous infrared galaxies (ULIRGs), which have column 
densities of Nh ~ 10 23 — 10 24 cm~ 2 . We now test this 
idea. In Figure 2, we plot the source-frame column den- 
sity probability distributions of the four DOAs in Table 
1 of [5] (see below). Redshifts are not known for three 
of these bursts, so we adopt z = 1 (top panel), 2 (middle 
panel), and 3 (bottom panel) for these three bursts. Given 
the redshifts that have been measured for bursts to date, 
the probabilities that all three of these bursts have z > 1, 
2, and 3 are 0.27, 0.013, and 0.0016, respectively. Even 
in the very unlikely event that all three of these bursts 
have z ~ 3, these four probability distributions are not 
consistent with the expected column density distribution 
for bursts that occur in the nuclear regions of ULIRGs, 
but are consistent with the expected column density dis- 
tribution for bursts that occur in molecular clouds [5]. 

Another possibility is that burst host galaxies are very 
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FIGURE 3. R-band magnitude or limiting magnitude, scaled 
to 18 hours after the burst, vs. 2 - 10 keV flux or limiting flux, 
scaled to 9 hours after the burst. The solid ellipse is from the 
best-fit model of §5 of [10] to all of the data, and the dotted 
ellipse is from the best-fit model to all of the data excluding 
GRB 970828. The ellipses should encompass sa 90% of the 
data. Clearly, GRB 970828 does not dominate the fit, and most 
of the DOAs are X-ray faint, but given these X-ray fluxes, the 
DOAs are optically fainter than the expectation from changes 
in collimation angle and distance alone (solid line). From [10]. 



dusty. However, if this were the case, one would ex- 
pect the X-ray bright afterglows to be DOAs as often as 
the X-ray dim afterglows, which is not the case (Figure 
3). However, Figure 3 is naturally explained by circum- 
burst extinction: The X-ray bright afterglows are proba- 
bly strongly collimated, and consequently the burst and 
optical flash are more likely to burn through the optical 
depth of the circumburst cloud; the X-ray dim afterglows 
are probably weakly collimated, and consequently the 
burst and optical flash are less likely to do so [10]. 

High redshift effects: Lyman limit absorption in the 
source frame, absorption by the Lya forest, absorp- 
tion by excited molecular hydrogen in the circumburst 
medium [14], and source-frame extinction by the FUV 
component of the extinction curve (e.g., [3]) could all 
result in DOAs if at sufficiently high redshifts. However, 
[12] show that unless the burst history of the universe dif- 
fers dramatically from the star-formation history of the 
universe, the redshift distribution of the bursts should be 
fairly narrowly peaked around z w 2, primarily because 
there is very little volume in the universe at low and high 
redshifts. [15] model the detection efficiency functions 
of BeppoSAX and IPN, and show that these satellites 
should detect even fewer bursts at high redshifts, push- 



ing the expected typical redshift for bursts detected by 
these satellites down to the observed value of about one. 
However, based on their variability redshift estimates of 
220 BATSE bursts, [16] (see also [17]) find that the burst 
history of the universe might differ dramatically from the 
star-formation history of the universe, with very many 
more bursts at high redshifts. However, using similar 
variability redshift estimates [4] for 907 BATSE bursts, 
[18] find that only w 15% of bursts above BATSE's de- 
tection threshold have z > 5 (and if the luminosity func- 
tion is evolving, far fewer bursts below BATSE's detec- 
tion threshold have z > 5). Consequently, the above high 
redshift effects probably affect <; 10% of the bursts in our 
BeppoSAX- and IPN-dominated sample. Furthermore, 
[13] find that the variability redshift estimates for all of 
the DOAs for which high resolution BATSE light curves 
are available have z < 5. 

A number of ad hoc hypotheses are also addressed in 
[10]. 



Evidence for a Large, Massive Cloud Origin 
of DOAs 

Neglecting for a moment extinction exterior to the 
circumburst cloud, we now consider the optical depth to a 
burst that is embedded a distance r within its circumburst 
cloud, along the line of sight. Prior to sublimation and 
fragmentation [7,8,9], the column to the burst is optically 
thin, which we take to mean x <; 0.3, if the hydrogen 
column density Nh 5 x 10 20 cm~ 2 , and optically thick, 
which we take to mean x J> 3, if Nh J> 5 x 10 21 cm~ 2 . 
We mark column densities with dotted lines in Figure 
4. Since sublimation and fragmentation burn through w 

lOL^ 2 pc of optical depth [7,8,9], where L = 10 49 t,49 erg 
s _1 is the 1 - 7.5 eV isotropic-equivalentpeak luminosity 
of the optical flash, the post-sublimation/fragmentation 
column to the burst is optically thin if r is less than this 
distance or Nh <, 5 x 10 20 cm~ 2 , and optically thick if 
r is greater than this distance and Nh ;> 5 x 10 21 cm~ 2 . 
We show this in Figure 4 by plotting x = 0.3 and 3 for 
L49 = 0.1 (thin curves, left), 1 (thick curves), and 10 (thin 
curves, right; see [7] for details). 

Consequently, modulo the value of L, bursts that oc- 
cur to the lower left of the x = 0.3 curve have either rel- 
atively unextinguished optical afterglows, or optical af- 
terglows that are extinguished by dust elsewhere in the 
host galaxy, or in our galaxy. Bursts that occur to the up- 
per right of the x = 3 curve have highly extinguished af- 
terglows. Since circumburst extinction appears to be re- 
sponsible for most of the DOAs, and Galactic extinction 
and host galaxy extinction unrelated to the circumburst 
medium account for no more than perhaps a few of the 
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FIGURE 4. Post-sublimation/fragmentation optical depth x 
to a burst that is embedded a distance r, along the line of sight, 
within a cloud of constant hydrogen density The three pairs 
of solid curves mark x = 0.3 (lower left) and 3 (upper right) for 
L49 =0.1 (left), 1 (center), and 10 (right). For a given L, bursts 
that occur to the lower left of the x = 0.3 curve have afterglows 
that are relatively unextinguished by the circumburst cloud, and 
bursts that occur to the upper right of the x = 3 curve have 
highly extinguished afterglows. The dotted lines mark constant 
hydrogen column densities Nh- From [7]. 

DOAs detected to date, we find that most DOAs have 
r £ 10L% 2 pc and N H ^5x 10 21 cm" 2 . 

We now estimate the sizes and masses of the circum- 
burst clouds of the DOAs. Taking the clouds to be spher- 
ical, and taking the bursts to be located at the cloud cen- 
ters, both of which are reasonable approximations on av- 
erage, we find that most DOAs occur in clouds of radius 

R > 10L% 2 pc and mass M £ 3 x 10 5 L 49 M (Figure 4 
of [7]). Clouds of this size and mass are typical of giant 
molecular clouds (e.g., [19]), and are active regions of 
star formation. 



EVIDENCE FOR A MOLECULAR 
CLOUD ORIGIN OF BURSTS 

In this section, we first consider the equilibrium prop- 
erties of clouds, and show that clouds of these sizes 
and masses, and hence the circumburst clouds of DOAs, 
cannot be diffuse, but rather are probably molecular, in 
which case they are active regions of star formation. 
Next, we show that Galactic-like molecular clouds have 
sufficiently high column densities to make bursts dark 



at optical wavelengths, and that sublimation of dust by 
the optical flash and fragmentation of dust by the burst 
and afterglow can leave a comparable number of bursts 
- bursts that are sufficiently near the earth-facing side of 
the cloud and/or in sufficiently small clouds - relatively 
unextinguished at optical wavelengths, which is in agree- 
ment with what is observed. Finally, we model and con- 
strain the column density distribution of the bursts with 
detected optical afterglows, and show that these bursts 
are also consistent with a molecular cloud origin, and are 
not consistent with a diffuse cloud origin. Consequently, 
we find that all but perhaps a few bursts, dark or other- 
wise, probably occur in molecular clouds. 



Evidence for a Molecular Cloud Origin of 
DOAs: Cloud Equilibrium, Gravitational 
Collapse, and Fragmentation-Driven Bursts 
of Star Formation 

We first consider the equilibrium properties of diffuse 
clouds: i.e., gas that is bound by pressure equilibrium 
with a warm or hot phase of the interstellar medium. 
Specifically, [20] considers the equilibrium properties of 
a uniformly magnetized, isothermal, non-rotating diffuse 
cloud, and introduces factors c\ = 0.53 and C2 = 0.60 that 
correct for the fact that in equilibrium the cloud will not 
be of uniform hydrogen density n#, but will be centrally 
condensed, and will not be spherical with radius R, but 
will be moderately flattened along the direction of the 
magnetic field. [20] shows that clouds with R > R m are 
unstable to gravitational collapse, where 
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M = 47tp/?^/3 is the mass of the cloud, p = /jf/m//n//, 
IAh = 1-87 is the mean molecular weight per hydrogen 
atom, 

3.15c 2 (M7/;) 4 



G 3 M 2 [l-(M t /M)2/3]3 



(4) 



T is the temperature of the cloud, /j = 1 .44 is the mean 
molecular weight, 



0.0236(ciB) 3 

G 3 /2p2 



(5) 



and B is the strength of the magnetic field within the 
cloud. We plot R,„(n H ) for T = 10 and 100 K and B = 3 
juG in Figure 5, and for T = 30 K and B = 1 and 10 /jG in 
Figure 2 of [5]: Most Galactic clouds have 10 < T < 100 
K, and B = 3 /jG is typical of the Galactic interstellar 
medium. 




FIGURE 5. Post-sublimation/fragmentation optical depth x 
to a burst that is embedded a distance r, along the line of sight, 
within a cloud of constant hydrogen density ng, along the line 
of sight (solid curves, from Figure 4), and diffuse cloud stability 
criteria (dashed curves). The dashed curves mark the maximum 
radius for which a diffuse cloud is stable against gravitational 
collapse, for cloud temperatures T = 10 and 100 K and cloud 
magnetic field strength B = 3 /jG. From [5]. 



Also in Figures 5 and 2 of [5], we plot the curves 
of constant post-sublimation/fragmentation optical depth 
X from Figure 4. Above, we argue that most of the 
DOAs occur to the upper right of the x = 3 curve, 
which implies relatively large sizes and high densities 
for their circumburst clouds. We show in Figures 5 and 
2 of [5] that diffuse clouds of these sizes and den- 
sities are unlikely to be in equilibrium, or even near 
equilibrium, for typical temperatures and magnetic field 
strengths. Such clouds would be unstable to gravita- 
tional collapse, resulting in the collapse and fragmenta- 
tion of the cloud until star formation re-establishes pres- 
sure equilibrium within the fragments, and the fragments 
are bound by self-gravity: This is precisely the structure 
of molecular clouds (e.g., [19]). Furthermore, molecu- 
lar clouds regularly have sizes R ;> 10 pc and masses 
M^3x 10 5 M , and are active regions of star forma- 
tion. Consequently, DOAs probably occur in collapsing 
diffuse clouds/forming molecular clouds, and/or molecu- 
lar clouds that re-established equilibrium some time ago. 
We confirm that molecular clouds have sufficiently high 
column densities to make a large fraction of bursts dark 
at optical wavelengths below. 

[19] find the mass scale of the molecular cloud frag- 
ments to be on the order of a few solar masses, which 



supports the idea that they are supported against further 
collapse and fragmentation by star formation. Since the 
lifetimes of massive stars are comparable to the free fall 
time of the progenitor cloud for these cloud masses, the 
first generation of massive star formation and supernovae 
should occur on this short timescale: <~ 10 7 — 10 8 yr. 
Consequently, DOAs are probably a byproduct of this 
burst of star formation if the molecular cloud formed re- 
cently, and/or the result of lingering or latter generation 
star formation if the molecular cloud formed some time 
ago. 

Consistency Check: The Column Density 
and Optical Depth Distributions of Bursts 
in Galactic-Like Molecular Clouds 

We now consider the column density and optical depth 
properties of bursts in Galactic-like molecular clouds. 
First, we compute the mean column densities of the 273 
molecular clouds in the sample of [19], which was con- 
structed using data from the Massachusetts-Stony Brook 
CO Galactic Plane Survey. The mean column density of 
a cloud is given by Nh = Mj (/j//m//A), where M is the 
virial mass ([19] show that the clouds are in or near virial 
equilibrium), A = 1 1 . 56 (D tan y/ a iOb) 2 is the effective 
cross section of the cloud [19], D is the distance to the 
cloud, and O/ and Ob are the angular sizes of the cloud. 
We plot the distribution (dotted histogram) of mean col- 
umn densities in Figure 6, and confirm the finding of [19] 
that hhNh is narrowly peaked around w 170 M pc~ 2 , 
independent of cloud mass and size. 

Next, we correct this distribution for a number of ef- 
fects. First, since bursts more likely trace cloud mass than 
cloud number, we weight these mean column densities 
by cloud mass, and renormalize the distribution (dashed 
histogram). Since molecular cloud column densities ap- 
pear to be fairly independent of cloud mass, the dis- 
tribution is not significantly changed. Also, we weight 
the M < 1 x 10 4 M© clouds by an additional factor of 
(M/7 x 10 4 M©) -3 / 2 to correct for an undercounting of 
low mass clouds on the far side of the Galaxy [19]. How- 
ever, this affects the mass-weighted mean column density 
distribution negligibly. 

Finally, we correct for geometrical effects: (1) bursts 
occur in their circumburst clouds, not behind them; and 
(2) molecular clouds are centrally condensed. First, we 
adopt the same density distribution for the clouds that 
[19] adopted to compute the virial masses: p °< r~ x 
within the effective radius of the cloud, and p = be- 
yond this radius. This density distribution results in a sur- 
face brightness profile that is similar to what is observed 
[19]. Next, we trace this density distribution with 10 4 
randomly-placed bursts for each cloud, and compute for 
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FIGURE 6. Expected column density distribution n(logA'//) 
(log-normalized to one) for bursts that occur in Galactic- 
like molecular clouds (solid histogram). The dotted histogram 
marks the mean column density distribution of the Galactic 
molecular clouds of [19], which we compute from their mass 
and size measurements. The dashed histogram marks the mass- 
weighted mean column density distribution, which we also cor- 
rect for an undercounting of low mass clouds on the far side of 
the Galaxy (a negligible correction). The solid histogram marks 
this distribution after correcting for the fact that bursts occur in 
their circumburst clouds, not behind them, and for the fact that 
molecular clouds are centrally condensed. From [5]. 

each burst the distance and column density to the earth- 
facing side of the cloud. The effect of placing bursts in 
the clouds, as opposed to behind the clouds, is to lower 
the mean column densities on average by a factor of two, 
and sometimes (when bursts occur near the earth-facing 
side of the cloud) considerably more. However, this is 
partially offset by the central condensation of the clouds, 
which places more bursts in above average density (and 
column density) environments. We plot the final distribu- 
tion - the expected column density distribution for bursts 
in Galactic-like molecular clouds - also in Figure 6 (solid 
histogram). 

Lastly, we confirm that Galactic-like molecular 
clouds, the column densities of which appear to be fairly 
universal, at least within the Galaxy, have sufficiently 
high column densities to make a large fraction of bursts 
dark at optical wavelengths. To this end, we plot in 
Figure 7 the expected distribution (1 and 2 o dashed 
contours) of column densities and distances to the 
earth-facing side of the cloud for bursts in Galactic- 
like molecular clouds. Also in Figure 7, we plot the 
curves of constant post-sublimation/fragmentation op- 
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FIGURE 7. Expected distribution of column densities and 
distances to the earth-facing side of the cloud for bursts that 
occur in Galactic-like molecular clouds (dashed contours), and 
the curves of constant post-sublimation/fragmentation optical 
depth from Figures 4 and 5 (solid curves). The thick dashed 
contour contains 68.3% of the distribution, and the thin dashed 
contour contains 95.4% of the distribution. The dotted lines 
mark constant hydrogen column densities. From [5], 



tical depth from Figures 4 and 5, and lines of constant 
column density. It appears that for low values of the 
isotropic-equivalent peak luminosity of the optical flash 
(L49 <; 0.1), Galactic-like molecular clouds can make 
a large fraction of bursts dark at optical wavelengths. 
However, even for very low values of L, a reasonable 
fraction of bursts would be left relatively unextinguished 
at optical wavelengths. Since comparable fractions of 
DOAs and detected optical afterglows are observed 
(about 2/3 and 1/3, respectively; e.g., [1,2]), this sug- 
gests that the bursts with detected optical afterglows 
might also occur in molecular clouds. We confirm this 
by modeling and constraining the distribution of col- 
umn densities, measured from absorption of the X-ray 
afterglow, of the bursts with detected optical afterglows 
below. 



Evidence for a Molecular Cloud Origin for 
Bursts with Detected Optical Afterglows: 
The Column Density Distribution 



pected distribution for bursts in Galactic-like molecu- 
lar clouds (solid histogram of Figure 6). To this end, 
we have reviewed the literature, and list in Table 1 
of [5] observer-frame column densities and uncertain- 
ties that have been measured from the X-ray after- 
glows of 15 bursts, 11 of which have detected op- 
tical afterglows (this information is not yet available 
for the vast majority of the DOAs). Before modeling 
these data, we convert each measurement Nhj and up- 
per and lower 1-a uncertainty G„.; and O/ , into proba- 
bility distributions, first in the observer frame, and then 
in the source frame by correcting the observer-frame 
probability distribution for the Galactic column density 
along the line of sight, and the redshift of the burst 
if known. If a u j = 0/ ]( (or a/ = Nhj), we take the 
observer-frame probability distribution to be a Gaussian 
of mean Nn,i and standard deviation a u /. Pi{NH, bs) — 
G{NHpbs,NH,i,<3u,i)- Otherwise, we take the observer- 
frame probability distribution to be a skewed Gaussian: 

Pi(N H ,obs) = G{Nl lobs ,Nl li ,o s i )dN s Hob JdN H .ob s , where 
s is the skewness parameter, and s and o, are given by 
solving 



= (N H ,i+o u .iY 



= Nh 



-N° H 



(6) 



i.e., we take the observer-frame probability distribution 
to be a Gaussian in b , which reduces to a Gaus- 
sian in Nn,obs when a u j = a; )( (s = 1). The source-frame 
probability distribution pi{Nh) is then given by substi- 
tuting Nn^obs = Nh.mw + N H (l+zy 2 - 6 , where N H 
the Galactic column density along the line of sight (inter- 
polated 2 from the maps of [21]), Nh is the source-frame 
column density, and (1 +z)~ 2 ' 6 scales Nh to the observer 
frame (e.g., [22]). If multiple column density measure- 
ments are available for a burst, we average their source- 
frame probability distributions if the measurements were 
made from the same data, and we take the product of 
these distributions if the measurements were made from 
independent data (GRB 000926). We plot the source- 
frame probability distributions, peak-normalized to one, 
for the bursts with detected optical afterglows in the top 
panel of Figure 8. The dotted curves mark bursts of un- 
known redshift, in which case we fix z = 0: These dis- 
tributions can be thought of as fuzzy lower limits (see 
below). 

We model the column density distribution of the 
bursts with detected optical afterglows with a broken 
power law. We use a three-parameter function for flex- 
ibility, and a simple function (as opposed to, e.g., the 
more elegant skewed Gaussian) so the inner integral of 
Equation (8) can be evaluated analytically (otherwise, 



In this section, we model and constrain the column 
density distribution of the bursts with detected opti- 
cal afterglows, and compare this distribution to the ex- 



2 We use the n# FTOOL at http://heasarc.gsfc.nasa.gov/cgi- 
bin/Tools/w3nh/w3nh.pl. 



1.2 




N H [cm" 2 ] 



distribution (solid curve), which is also log-normalized, 
and the l-o uncertainty in this distribution (dotted 
curves) in the bottom panel of Figure 8. We find that the 
column density distribution of the bursts with detected 
optical afterglows peaks around 10 22 cm~ 2 , and spans 
a factor of a few to an order of magnitude or so. These 
findings verify and improve upon the earlier finding of 
[6] that these column densities span 10 22 - 10 23 cm . 

Lastly, we compare the fitted column density distri- 
bution of the bursts with detected optical afterglows to 
the expected column density distribution for bursts in 
Galactic-like molecular clouds (solid histogram, replot- 
ted from Figure 6), and find these distributions to be con- 
sistent within the l-o uncertainties. Furthermore, the fit- 
ted distribution is not consistent with the expectation for 
bursts that occur in diffuse clouds, the canonical column 
density of which is a few times 10 20 cm~ 2 . Consequently, 
we find that all but perhaps a few bursts, dark and other- 
wise, probably occur in molecular clouds. 



FIGURE 8. Top panel: Source-frame column density proba- 
bility distributions pi(Nh) (peak-normalized to one) for the 11 
bursts with detected optical afterglows in Table 1 of [5]. The 
dotted curves mark bursts of unknown redshift, for which we 
fix z = 0: These distributions can be thought of as fuzzy lower 
limits. Bottom panel: Likelihood-weighted average model dis- 
tribution n(logNH) (solid curve), and l-o uncertainty in this 
distribution (dotted curves). The solid histogram is the expected 
column density distribution for bursts that occur in Galactic- 
like molecular clouds from Figure 6. These distributions are 
consistent within the l-a uncertainties. From [5]. 



the error analysis would be computationally taxing). 
We note that the broken power law is log-normalized 
[J^n^ogN^dlogNn = 1] to facilitate comparison of 
the fitted distribution to the solid histogram of Figure 6 
(see below). 

We fit the model «(logAf//) to the data pi(Nh) using 
Bayesian inference (e.g., [3]). The likelihood function is 
given by Equation (1), where N = 11 is the number of 
bursts with detected optical afterglows in our sample, and 



Li 



P„Pi{N H )n{N H )dN H 
I^Pi(N' H ) \On(N H )dN H 



dN' H 



(Zi known) 
(Zi unknown) 
(7) 



(e.g., [10]), or equivalently 



I-oo Pi (NH)n(logN H )d \ogN H 
!oPi( N H) [flogN' H n(\ogN H )dlogN H 



dN' H 



since n(logA r //) 
log(a/cm~ 2 ) = 



for N H < 0. We find 



22.16 



-0.20 
-0.33' 



arctanb 



0.92 



(8) 
that 

-0.38 
-0.35 



(b w 1.31^|), and arctanc = -1.46+J^ [c ^ -1.33 
(1 a)]. We plot the likelihood-weighted average model 
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